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The paper derives closed form expressions for transient and steady-state 
operation of a DC-DC converter with single switched capacitor. To this end, 
the result of each switching is considered as a point in the iterative process, 
and the function between the points is reconstructed. As opposed to the 
commonly accepted approach, when each of the topologies is approximated 
by a first order circuit, the proposed analysis is carried out for second order 
circuits. This allows obtaining the waveform of output voltage ripple and 
paves the way to more accurate calculation of equivalent resistance. The 
obtained analytical expressions were verified by simulations and an excellent 
agreement between the results was found. 
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1. INTRODUCTION 

Switched capacitor converters (SCCs) are favored in some applications due to low EMI and 
compatibility with IC technology. Over the past few decades, ongoing research has shifted towards 
sophisticated SCCs with large number of capacitors and advanced control circuits. However, only a few 
studies use analytical methods. The analysis presented in this paper is based on the method of difference 
equations [l]-[3] which allows finding the solution for the transient and steady-state operation. Although this 
method is unified, it was applied previously only to the switched inductor converters [4]-[6]. However, some 
analytical methods for analysis of switched capacitor circuits that bear resemblance of the proposed one, 
were presented in [7]-[10]. 

Let us consider the SCC shown schematically in Figure 1(a). It comprises four switches S ± ... S 4 with 
on-resistances r t ...r 4 . The corresponding pair of switches is turned on/off by two non-overlapping clocks (p 1 
and cp 2 shown in Figure 1(b). Thus, during t x the capacitor C ± is charged by V in through S l9 S 3 and then, 
during t 2 , is discharged to the load through S 2 , S 4 . Thus, we have two topologies, which are considered 
separately in Sections II. Applying to each topology the Kirchhoff s voltage and current laws (KVL and 
KCL), we write a system of two first-order differential equations. These equations are then solved using the 
Laplace transform. The solution for the first topology defines the initial conditions for the second one. This 
enables us to compose a system of two first-order difference equations, which is solved using the Z- 
transform. Thus, for a given period, n, we know the initial values of state variables and functions according 
to which these variables change. This is the solution is closed form. 
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Figure 1. (a) Considered SCC (b) and two non-overlapping clocks cp t and cp 2 


2. DIFFERENTIAL EQUATIONS 
2.1. First topology 

The switches S ± and S 3 is shown in Figure 1(a) are turned on during nT < t <nT + t l9 where n is 
the number of periods. Thus, in the first topology as shown in Figure 2, R ± = r t + r 3 . 
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Figure 2. First topology of the considered SCC 
The KVF and KCF equations for this circuit are 
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Fet us write (1) in the Faplace domain using 
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such that 
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where 
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The inverse Laplace transform of (4) is: 
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Using t = t ± and 

we can rewrite (7) as: 

Vaiti) = a + b(V :r , - V C0 (n}) + cV ci in) 
VcJt;) = d(V in - V C1 (n)) -|- eV Co (;n) 


where 



2.2. Second topology 

The switches S 2 and S 4 in Figure 1(a) are turned on during nT + t 4 < t < (n + 1 )T, 
second topology as shown in Figure 3 R 2 = r 2 + r 4 . 



Figure 3. Second topology of the considered SCC 
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such that in the 
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The KVL and KCL equations for this circuit are: 
Fci = — + f/ CG 
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l-h = h 
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dV cl 
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Using (2), we write (11) in the Laplace domain: 

| + 2}E lC i - V Co = ■RjQP r C j(n) 
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The solution of (12) is: 
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The inverse Laplace transform of (13) is: 
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Since for the second topology t t < t < t 2 , the initial conditions will be L cl (n) — Lci(U) and 
VcoW = Vco(.h ), whereas t = t 2 - t ± . 

3. DIFFERENCE EQUATIONS 

Substituting (14) and (15) into (16), we obtain the recurrent equations: 


where 
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Now, substituting (9) into (17), we obtain the difference equations: 

(V cl {n + 1) - A + BV Ca (n) - GV ci (n) 

\V co in +1)=E+ FV C0 (n) + HV cl (n) 

where 

A=(fb + gd)V in -fa; B = ge - fb; G = fc - gd 
E = (hb- f- J«OP in -I- ha; F = ke - hb; H = he - kd 

Before we proceed to the solution of (19), let us consider the steady-state operation, where 

Vain + 1} = Vci(n>' V Co (n -f l) = V Co (n) 

Substituting (21) into (19), we have 

f(c - i)F C3 («) + BV Co (n) = -A 
WaW + tF- D^OO = -S 

The solution of (22) is: 
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Vcl ^ ~ (f - l)(G - 1) - BH' Vca< " n> ~ (F- 1)(G - 1) - BH 

Let us write (19) in the Z-domain using 

Z{K) = A F —and Z{f(n + l)} = zF(n) - f(0) 
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Since Error! Reference source not found, is represented in the form of partial fractions, we can apply the 
inverse Z-transform to each 
term separately: 
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4. SIMULATION RESULTS 

To verify the obtained analytical expressions (9), (17) and (28), the circuit shown in Figure 4 was 
simulated in PSIM 9.0. Since the PSIM bidirectional switches have zero on-resistance, two external resistors 
corresponding to R ± and R 2 in Figure 2 and Figure 3 were added. The parameters of the circuit in Figure 4 
are as follows: V in = 10F, F C1 (0) - Fco(O) - OF, h = t 2 = 5/rs, R ± = R 2 = 1H, C ± = 10 p.F, R 0 = 100H 
and C Q = 100/rF. 



Figure 4. Simulation circuit for the voltage-halving SCC 


Since R ± = R 2 and t ± = t 2 , the expressions (10) and (18) are reduced to: 
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These constants are substituted into (20) and then into (28), which sets the initial conditions for (9) and (17). 
The voltages across C ± and C 0 during the first ten periods (0 < t < 107) is shown Figure 5, which compares 
the MathCAD calculation and PSIM simulation. 
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(a) 



(b) 


Figure 5. Voltages across C ± and C 0 for 0 < t < 107, MathCAD (a) and PSIM (b) The horizontal scale is 

10 fis/div, i.e. each division is 7 


The voltages in Figure 5 were measured at the points nT and are given in Table 1 along with the 
relative error, £. 


Table 1. Measured values of the voltages in Figure 5 


n 

Vci 

Calc. 

(n) 

Simul. 

£[%] 

Vco 

Calc. 

(n) 

Simul. 

£[%] 

1 

2.5146 

2.5046 

0.398 

0.5174 

0.5119 

1.074 

2 

3.4986 

3.4937 

0.140 

0.9542 

0.9449 

0.984 

3 

3.9147 

3.9144 

0.008 

1.3381 

1.3259 

0.920 

4 

4.1172 

4.1193 

0.051 

1.6815 

1.6670 

0.870 

5 

4.2366 

4.2397 

0.073 

1.9909 

1.9747 

0.820 

6 

4.3211 

4.3254 

0.010 

2.2706 

2.2531 

0.777 

7 

4.3890 

4.3943 

0.121 

2.5236 

2.5054 

0.726 

8 

4.4474 

4.4534 

0.135 

2.7526 

2.7339 

0.684 

9 

4.4991 

4.5055 

0.142 

2.9601 

2.9410 

0.649 

10 

4.5454 

4.5455 

0.002 

7.1478 

7.1282 

0.627 


Figure 6 compares the MathCAD calculation and PSIM simulation for 0 < t < 1007. Note that at 
t ^ 607 the SCC reaches the steady-state 
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Figure 6. Voltages across C x and C 0 for 0 < t < 1007, MathCAD (a) and PSIM (b). 
The horizontal scale is 100 fis/div, i.e. each division is 107 


The steady-state voltages are shown Figure 7. Their discrete values were calculated by (23) and then 
substituted as the initial conditions into (9) and (17). 



Figure 7. Steady-state voltages across C ± and C 0 for 957 < t < 1007 
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5. CONCLUSIONS 

Based on the method of difference equations the closed form expressions for the voltages across the 
capacitors in the voltage-halving SCC were derived. The solution of these equations allows us to predict the 
SCC behavior in both the transient and steady-state operation. That is for a given period, n, we know the 
initial values of the voltages and functions according to which these voltages change. The obtained 
expressions were verified by simulations. As evident from Table , the deviation between the theoretical and 
simulation results does not exceed 1.1%. The used method however, is very complex even in the case of the 
two-phase SCC and its extension to the multi-phase SCC will apparently require some special assumptions. 
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